Orbital and superorbital variability and their coupling in X-ray binaries 



Andrzej A. Zdziarski, Juri Poutanen, Askar Ibragimov, Marek Gierlinski, '° and Linqing Wen' 



^6,7,8 



00 

o 

O 
(N 

D 
(N 



o 



o 

G\ 
O 
oo 
O 



^Centrum Astronomiczne im. M. Kopernika, Bartycka 18, 00-716 Warszawa, Poland 
^Astronomy Division, Department of Physical Sciences, PO Box 3000, FIN-90014 University of Oulu, Finland 
^Kazan State University, Astronomy Department, Kremlyovskaya 18, 420008 Kazan, Russia 
'^Department of Physics, University of Durham, Durham DHl 3LE, UK 
^Obserwatorium Astronomiczne Uniwersytetu JagieUohskiego, Orla 171, 30-244 Krakow, Poland 
^Max-Planck-Institut fiir Gravitationsphysik, Albert-Einstein-lnstitut, Am Miihlenberg 1, D-14476 Golm, Germany 
^Division of Physics, Mathematics and Astronomy, Caltech, Pasadena, CA 91125, USA 
^School of Physics, University of Western Austraha, Crawley, WA 6009, Australia 

E-mail (AAZ): aaz@camk.edu.pl 

Abstract 

We review X-ray flux modulation from X-ray binaries on time scales corresponding to the orbital period 
and those at longer time scales (so called superorbital). Those modulations provide a powerful tool to 
constrain geometry of the accretion flow. The most common cause of the superorbital variability appears 
to be precession. We then discuss two specific examples of discoveries of a coupling between the two types 
of variability and their physical interpretation. One is Cyg X-1, a black-hole binary with a high-mass 
companion, in which case we find the presence of an accretion bulge formed by collision of the stellar wind 
with the outer edge of the precessing accretion disc. The other is 4U 1820-303, a neutron star accreting 
from a low-mass white dwarf, in which case we interpret the superorbital variability as accretion rate 
modulation induced by interactions in a triple stellar system. Then, the varying accretion rate leads to 
changes of the size of the accretion bulge in that system, obscuring the centrally-emitted X-rays. 

Key WORDS: accretion, accretion discs — stars: individual: (Cyg X-1, HDE 226868, 4U 1820-30) — 
X-rays: binaries — X-rays: stars 



1. Introduction 

A large number of X-ray binaries show X-ray flux period- 
icities at their orbital period. The current comprehensive 
list of such objects is given by Wen et al. (2006). Using 
strict criteria for the statistical significance of the pres- 
ence of a modulation, they find 33 secure orbital periods 
in the RXTE / ASM data. Due to the fimitation of the 
ASM sampling, they were unable to detect the 685-s X- 
ray binary period of the ultracompact low-mass X-ray 
binary (LMXB) 4U 1820-303, which, however, is clearly 
detected, e.g., in RXTE /VGA data (e.g., Zdziarski et 
al. 2007b, hereafter Z07, and references therein). Even 
faster X-ray modulations, at 569 s and 321 s, are found in 
ROSAT data for RX J1914.4-I-2456 (Cropper et al. 1998) 
and RX J0806.3-t-1527 (Israel et al. 1999), respectively. 
These two systems contain pairs of white dwarfs orbiting 
each other at the above respective period (Cropper et al. 
1998; Israel et al. 2002; Ramsay et al. 2002). 

In addition, a number of X-ray binaries show X-ray 
modulations at periods longer than their orbital peri- 
ods, so-called superorbital periodicity. Sood et al. (2007) 
state that 35 such systems have been reported in litera- 



ture. However, many of those claims remain uncertain, 
and Sood et al. (2007) selected 19 of them as secure. 
However, even some of those appear, in our opinion, un- 
certain. For example, they list the high-mass X-ray bi- 
nary (HMXB) Gen X-3 as having a superorbital period of 
140 d. On the other hand, the same object is presented 
by Raichur & Paul (2008) as a case of the lack of stable 
superorbital modulation. Another example is 4U 1916- 
053, which is listed as having a 199-d period in Sood et al. 
(2007) whereas Homer et al. (2001) (and later Wen et al. 
2006) find no trace of such periodicity. Wen et al. (2006), 
using a uniform analysis method, find 6 objects as having 
coherent superorbital modulations in the RXTE / ASM 
data, and further 5 with quasi-periodicities. Using their 
method with no pre-selection of data, they did not find a 
superorbital periodicity in the black-hole HMXB Cyg X- 
1. However, its '^150-d period manifests itself mostly in 
the hard spectral state, where its statistical significance 
is very high (e.g., Brocksopp et al. 1999a; Karitskaya 
et al. 2001; Ozdemir & Demircan 2001; Lachowicz et 
al. 2006; Ibragimov et al. 2007; Poutanen et al. 2008, 
hereafter PZI08). Another black-hole binary in which 



the superorbital modulation manifests itself only in the 

hard state appears to be the LMXB GX 339 4. in which 
it shows a ~220-d periodicity, whereas that modulation 
is absent in the soft state (Zdziarski et al. 2004). (See, 
e.g., Zdziarski & Gicrlihski 2004; Done et al. 2007 for 
reviews and physical interpretation of the spectral states 
of black- hole binaries.) 

The MAXI, an X-ray all-sky monitor mission on the 
Japanese Experiment Module of the International Space 
Station, is scheduled to be launched in 2009 March. Due 
to its unprecedented sensitivity of a few mCrab in a day 
covering most of the sky, the MAXI will monitor vari- 
ability of a large number of X-ray sources at much lower 
flux levels than is now possible with the RXTE/ ASM. It 
will thus be able to greatly enhance our knowledge of the 
variability of X-ray binaries, in particular of their orbital 
and superorbital periodicities. 

2. Physical causes of orbital and superorbital flux period- 
icities 

An orbital periodicity may be caused by a number of ef- 
fects. Here, we review also some effects affecting other 
wavelengths than X-rays. First, a source associated with 
a compact object in a binary (usually a HMXB) may be 
eclipsed by its companion, see a list in Wen et al. (2006). 
Second, a flux modulation may be caused by an optically- 
thick disc rim (which is highest at the point of impact of 
the gas stream from the inner Lagrangian point in case of 
a donor filling its Roche lobe), obscuring the disc and/or 
its corona (e.g.. White & Swank 1982; Hellier & Mason 
1989; Heinz & Nowak 2001), usually in LMXBs. This ob- 
scuration may lead to strong partial eclipses in so-called 
X-ray dippers. More generally, the disc and any associ- 
ated structures may depart from its axial symmetry due 
to the influence of the companion, which may cause an 
orbital modulation. Third, wind from a high-mass com- 
panion may absorb/scatter the emission from the vicinity 
of the compact object, and the degree of absorption will 
depend on the orbital phase. In the case of Cyg X-1, 
both X-ray and radio emission are modulated by this ef- 
fect, which modulations were modeled by, e.g.. Wen et 
al. (1999) and Szostck & Zdziarski (2007), respectively. 
Fourth, phase-dependent absorption (via photon-photon 
pair production) of high-energy 7-rays may occur in a 
photon field axially asymmetric with respect to the com- 
pact object, especially that of the stellar photons (e.g., 
Bednarek 1997; Dubus 2006). A fifth effect of the com- 
panion is reflection or reprocessing of the emission from 
around the compact object on the surface of the com- 
panion facing the compact object (Basko 1978). This 
effects appears to be responsible for, e.g., the UV flux 
modulation from the X-ray binary 4U 1820-303 (Arons 
& King 1993; Anderson et al. 1997). Finally, the opti- 
cal/UV emission of the companion will be modulated if 



its shape departs from the spherical symmetry by par- 
tially or fully filling its Roche lobe, which effect is seen 
in Cyg X-1, e.g., Brocksopp et al. (1999b). 

Then, there will be an intrinsic dependence of the 
emitted flux on the orbital phase if the orbit is ellip- 
tical. This leads, e.g., to periodic outbursts around the 
periastron of Cir X-1 (e.g., Parkinson et al. 2003) and 
Be/X-ray binaries (see, e.g., Coe 2000; Neguerucla 2004 
for reviews) in X-rays, and sometimes, at other wave- 
lengths. Also, some orbital flux modulation may be due 
to the Doppler effect, which is in principle observable 
(Postnov & Shakura 1987), but has not yet been de- 
tected in a binary. (Obviously, the Doppler effect leads 
to widely observed shifts of spectral lines from binaries.) 

The observed superorbital variability appears in most 
cases compatible with being caused by accretion disc 
and/or jot precession, which either results in variable ob- 
scuration of emitted X-rays as in Her X-1 (Katz 1973), 
or changes the viewing angle of the presumed anisotropic 
emitter, as in SS 433 (Katz 1980) or Cyg X-1 (e.g., La- 
chowicz et al. 2006; Ibragimov et al. 2007; PZI08), or 
both. The precession may be induced via tidal forces by 
the companion (Larwood 1998), with additional forces 
exerted by irradiation of the disc by the central X-ray 
source (e.g., Ogilvie & Dubus 2001). The only certainly 
known exception, in which the superorbital periodicity 
is clearly caused by modulation of the accretion rate, 
M (and thus not by a changing viewing angle of the 
source), appears to be 4U 1820-303 (Zdziarski et al. 
2007a; Z07, and references therein). The M modulation 
can be achieved in a triple system via a secular modula- 
tion of the eccentricity of the inner system, proposed for 
this source by Chou & Grindlay (2001), and calculated 
in detail by Zdziarski et al. (2007a). Another system 
possibly of this type is the HMXB 2S 0114-f-650, where 
Farrell et al. (2008) have argued for changes of the M 
based on the spectral variability of the system. 

3. Coupling between orbital and superorbital modulations 
A number of binaries show both orbital and superorbital 
modulations. Those currently known are LMC X-4, 2S 
0114+650, SMC X-1, Her X-1, SS 433, 4U 1820 303 and 
Cyg X-1. An interesting issue then is whether there is 
any dependence of the parameters of the orbital mod- 
iflation on the superorbital phase (or. similarly, on an 
average of the flux level). The shape of the proflle of 
the orbital modulation in Her X-1 was found to depend 
on its superorbital phase (Scott & Leahy 1999), which 
appears to be due to the shadowing effect of the precess- 
ing accretion disc and scattering in its wind in that sys- 
tem. Analogous dependencies of the shape of the orbital 
modulation on the average flux level have been found 
(Raichur & Paul 2008) in LMC X-4, SMC X-1, Her X-1, 
as well as in Cen X-3 (which object, however, shows a 
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(a) A view along the orbital plane, =0 
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Fig, 1, (a) The superorbital phase diagram for the ASIVl 1,5-3 keV 
channel (the unit of F is count s~^). (b) Comparison of the 
characterization of the 1.5—3 keV orbital modulation dependence 
on the superorbital phase using different methods. The crosses 
show the intrinsic fractional rms of the orbital modulation. The 
solid histogram gives the orbital variability amplitude fitted by sum 
of three harmonics. The dashed histogram gives the corresponding 
rms for the fitting function. The solid curves in (a— b) show the 
dependencies for a theoretical outflow model fitted to the data 
by PZI08. (c) The mean hardness ratio (5-12 keV/1.5-3 keV), 
and (d) the distribution of X-ray dips in the hard state. A visible 
offset of the ephemeris with respect to 4> = is due to the use of 
an ephemeris based also on older data (Lachowicz et al. 2006). 



chaotic superorbital behavior). 

The Be/X-ray binary LS I +61°303 shows orbital vari- 
ability in the radio, X-ray and TeV emission, and a su- 
perorbital variability of the peak radio flux during an or- 
bit (Gregory et al. 1999; Gregory 2002). Gregory (2002) 
found a marked dependence of the phase of the peak of 
the orbital radio modulation on the superorbital phase. 
The presence of such a dependence may be due to in- 
teraction of the pulsar in that system with a variable 
circumstcllar Be decretion disc (Gregory 2002; Zdziarski 
et al. 2008), but details remain unknown. 

Then, strong coupling between the orbital and su- 
perorbital modulations is seen in the HMXB Cyg X-1 
(PZI08) and in the LMXB 4U 1820-303 (Z07), which 
consists of a white dwarf accreting onto a weakly- 
magnetized neutron star. We review these two objects 
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Fig. 2. A bulge at the outer edge of a precessing inclined disc. The 
material in the bulge absorbs some of the X-ray emission origi- 
nating close to the disc center. The orbital modulation due to 
the bulge is seen to strongly depend on the superorbital phase. 
(Orbital modulation due to the stellar wind itself is not illustrated 
here for clarity.) A view along the orbital plane: (a) the superor- 
bital phase 4> = 0, when the disc is seen closest to edge-on and 
the effect of the bulge is strongest; (b) the opposite case of the 
superorbital phase <E> = 0.5. (c) A view from the top. 



in detail below in Sections 4.-5. 

4. The dependence of the orbital modulation on the su- 
perorbital phase in Cyg X-1 

Such a dependence was searched for, studied, and 
interpreted theoretically by PZI08. They used the 
RXTE / ASM data, in which both orbital and superor- 
bital modulations are clearly seen, see Section 1. The 
physical cause of the orbital X-ray modulation is bound- 
free absorption in the stellar wind, whereas the widely 
accepted underlying cause of the superorbital, '--^ISO-d, 
modulation is accretion disc precession, see Section 2. As 
found by Ibragimov et al. (2007), the resulting superor- 
bital flux variability is most likely due to an anisotropy of 
the X-ray emission from the accretion flow (rather than 
due to absorption, obscuration or scattering). 

PZI08 found the fractional amplitude of the orbital 
modulation strongly depends on the superorbital phase 
in the hard state, see Fig. la, b. Namely, the orbital 
modulation is strongest when the average flux level is 
lowest. The effect is visible in all of the three ASM chan- 
nels, but it is strongest in the lowest-energy channel, 1.5- 




Fig. 3. (a) The RXTE /ASM light curve of 4U 1820-303. The short arrows shows the times of detected X-ray bursts, (b) The PCA light curve 
based on Galactic Bulge scans, (c) The light curve of the bolometric flux calculated from the PCA/HEXTE data. See Z07 for details. 
The ASM and PCA light curves have been extended up to MJD '~54700. The vertical lines show the minima (the phase $ = 0) of the 
superorbital cycle according to the ephemeris of Chou ii Grindlay (2001), which ephemeris is seen to become increasingly inaccurate in 
recent epochs. 



3 keV, where the orbital modulation amplitude is highest 
(which reflects the energy dependence of the bound-free 
cross section averaged over the cosmic abundances). 

A theoretical interpretation of this effect proposed by 
PZI08 is as follows. The cause of the part of the or- 
bital modulation dependent on the superorbital phase is 
a structure at the outer edge of the accretion disc on the 
side of the companion. Such a structure is expected the- 
oretically. Formation of the accretion disc in Cyg X-1 by 
a focused wind (Gies & Bolton 1986) leads, most likely, 
to a condensation of the wind matter near the disc outer 
edge on the side of the companion in the form of a bulge, 
similar to the disc bulge inferred to be present in low- 
mass X-ray binaries, e.g., White & Holt 1982; White 
& Swank 1982; Parmar & White 1988). On the other 
hand, the bulge can also be formed (see, e.g., Boroson et 
al. 2001) by a shock wave in the wind when it encoun- 
ters the gravity of the companion, the disc, or a wind 
from the disc, which is also likely to be present. In any 
case, when the fast, > 1000 km s~^, wind is stopped, the 
density increases dramatically. 

Such a bulge is shown in Fig. 2. We see that at the 
minimum of the superorbital cycle (i.e., at the superor- 



bital phase of $ ~ 0), when the disc is most edge-on 
(Fig. 2a), the bulge crosses the line of sight to the X-ray 
source (which is close to the black hole) around the supe- 
rior conjunction, i.e., at the orbital phase of ~ 0, and it 
is away from the line of sight at the inferior conjunction, 
i.e., at the orbital phase of ~ 0.5. On the other hand, 
at the maximum of the superorbital cycle ($ ^ 0.5), the 
disc is most face-on (Fig. 2b), and the bulge is away from 
the line of sight at any orbital phase. Thus, the bulge 
introduces additional absorption only around $ ~ 0.5, 
just as observed. 

This model was put into the form of a set of equations 
and then fitted to the data by PZI08. They considered a 
number of models for the anisotropy of the X-ray emis- 
sion, leading to the superorbital modulation, following 
Ibragimov et al. (2007). An exponential absorption pro- 
file of the bulge and either an isotropic or focused wind 
were assumed. In all models, the total optical depth of 
the bulge measured from its center was found to be 2±1. 
The bulge center was found to be displaced from the line 
connecting the centers of the stars by ~ 25°. The tilt 
of the precessing disc was found to be compatible with 
~ 10°, at the assumed inclination of 40°. Prograde pre- 
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Fig. 4. (a) The intrinsic colour-colour diagram, as given by the ra- 
tio of the energy fluxes in the photon energy ranges of (9.7—16 
l<eV)/(6.4-9.7 keV) and (4.0-6.4 keV)/(3. 0-4.0 keV) for the hard 
colour and soft colour, respectively. The boxes show a division into 
8 spectral substates. (b) The corresponding dependence of the 
hard colour on the bolometric flux. From Z07. 



cession was strongly favored by the fits. 

PZI08 also have shown that both the X-ray spectral 
hardness and the frequency of X-ray dips (caused by ab- 
sorption by wind blobs, Balucihska-Church et al. 2000) 
increase towards the superorbital phase of 0, see Fig. 
Ic, d, which is also explained by the above model. Fur- 
thermore, the coupling of the orbital and superorbital 
modulations leads to appearance of asymmetric beat fre- 
quencies in the power spectrum, explaining the finding 
of such asymmetry by Lachowicz et al. (2006). 

5. The dependence of the orbital modulation on the spec- 
tral state in 4U 1820-303 
Such a dependence was searched for, studied, and in- 
terpreted theoretically by Z07. Long-term light curves 
of the system are shown in Fig. 3. We see a very pro- 
nounced superorbital variability, with the period of ~170 
d. The strongest argument for the intrinsic, M-driven, 
character of this variability is the occurrence of X-ray 
bursts (which generally take place at low accretion rates) 
only during the deep minima of the light curve, see Fig. 
3a. As proposed by Chou & Grindlay (2001), the physi- 
cal cause of the M modulation is quasi-periodically vary- 
ing eccentricity, e, of the binary induced by a third star 
in the system, via the so-called Kozai effect (Kozai 1962). 
This was calculated in detail by Zdziarski et al. (2007a), 
who have modeled such e evolution and calculated the 
resulting M variability in a Roche-lobe overflow model. 

The orbital modulation of the system was discovered 
in X-rays (Stella et al. 1987). Given the very short orbital 
period, it is not detectable in the ASM data, and there- 
fore Z07 used PCA/HEXTE observations of 4U 1820- 
303. Their light curve is shown in Fig. 3c. We see their 
coverage of the superorbital variability is rather sparse, 
and it has proven difficult to study the orbital modula- 
tion directly as a function of the superorbital phase. On 
the other hand, the superorbital cycle is closely related 
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average (which are negligibly small). From Z07. 



to the spectral variability of the system, going between 
the so-called atoll (at low M) and banana (at high M) 
spectral states (Bloser et al. 2000). Therefore, Z07 have 
divided the pointed X-ray observations into 8 spectral 
substates, with the M increasing from the state 1 to 8, 
as shown in Fig. 4. 

Then the orbital modulation was fitted within each 
substate by a sinusoidal dependence. Fig. 5 shows the 
results. We see that the PCA count rate increases with 
the state number (as expected). The fractional modula- 
tion amplitude very significantly increases from the atoll 
(1-2) to banana (3-8) states. Also, the offset phase in- 
creases. Z07 interpreted the results in terms of the size 
and location of the bulge at the disc rim (partially ob- 
scuring the central source) changing with the variable 
accretion rate. Z07 also found that the fractional mod- 
ulation was independent of the photon energy. This ar- 
gues for the part of the bulge obscuring the X-rays being 
almost completely ionized. More details and discussion 
of other options are given in Z07. 

6. Discussion 

We have reviewed orbital and superorbital modulation 
of X-ray binaries, concentrating on the X-ray spectral 
band. The physical causes of the orbital modulation are 
relatively well understood, see Section 2. The most com- 
mon cause of the superorbital modulation appears be 
precession of accretion disc. However, details remain un- 
clear, in particular what drives the precession, whether 
it is prograde or retrograde, and how a viscous disc can 



coherently precess. On the other hand, some objects, es- 
pecially 4U 1820-303, show very strong evidence of the 
superorbital variability being due to periodic changes of 
M, and not due to precession. However, physical details 
also here remain unclear. One possibility is the effect 
of a third star in a hierarchical triple (Chou & Grindlay 
2001; Zdziarski et al. 2007a). 

Then, the properties of the orbital modulation may de- 
pend on the superorbital phase. We have discussed two 
specific examples. In Cyg X-1, the orbital modulation is 
strongest at the lowest observed flux, which appears to 
be due to the presence of an accretion bulge at the outer 
edge of the precessing disc. The bulge goes into the line 
of sight when the disc is most edge-on, which then cor- 
responds to the lowest observed flux. An opposite effect 
is observed in 4U 1820-303. Here, the orbital modula- 
tion is weakest at the lowest observed flux. The reason 
for that appears that the size of the obscuring bulge in- 
creases with the increasing M, periodically changing in 
that system (unlike Cyg X-1, where most likely there is 
no periodic M modulation). 
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